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Abstract 

In the search for bottom squark (sbottom) in SUSY at the LHC, the common practice has 
been to assume a 100% decay branching fraction for a given search channel. In realistic MSSM 
scenarios, there are often more than one significant decay modes to be present, which significantly 
weaken the current sbottom search limits at the LHC. On the other hand, the combination of the 
multiple decay modes offers alternative discovery channels for sbottom searches. In this paper, we 
present the sbottom decays in a few representative mass parameter scenarios. We then analyze 
the sbottom signal for the pair production in QCD with one sbottom decaying via b —)■ bxi, by^, 
and the other one decaying via b —)• txf ■ With the gaugino subsequent decaying to gauge bosons 
or a Higgs boson X2 ^Xii and xf ^^X?) we study the reach of those signals at the 
14 TeV LHC with 300 fb“^ integrated luminosity. For a left-handed bottom squark, we find that 
a mass up to 920 GeV can be discovered at 5a significance for 250 GeV < m^o < 350 GeV, or 
excluded up to 1050 GeV at the 95% confidence level for the h channel (/i > 0); similarly, it can 
be discovered up to 840 GeV, or excluded up to 900 GeV at the 95% confidence level for the Z 
channel (^ < 0). The top squark reach is close to that of the bottom squark. The sbottom and 
stop signals in the same SUSY parameter scenario are combined to obtain the optimal sensitivity, 
which is about 150 GeV better than the individual reach of the sbottom or stop. For a right- 
handed bottom squark with bb* —>• bxi, tx^ channel, we find that the sbottom mass up to 880 
GeV can be discovered at 5a significance, or excluded up to 1060 GeV at the 95% confidence 
level. 
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I. INTRODUCTION 


The outstanding performance of the Large Hadron Collider (LHC) at CERN had led 
to the milestone discovery of the Higgs boson predicted by the Standard Model (SM). 
High energy physics has thus entered a new era in understanding the nature of electroweak 
symmetry breaking. “Naturalness” argument for the Higgs boson mass implies new physics 
associated with the SM Higgs sector not far above the TeV scale [1-3]. The LHC Run-2 
with higher energy and higher luminosity will certainly extend the horizon to seek for new 
physics. Among the new physics scenarios, the weak scale Supersymmetry (SUSY) remains 
to be the most attractive option because of the accommodation for a light Higgs boson, 
the natural dark matter candidate, and the possibility for gauge coupling unification. In 
preparing to exploit a large amount of the incoming data from the LHC experiments, it is 
thus of priority to embrace the SUSY searches in a comprehensive way. 

While the top squark (stop t) sector might be the most relevant supersymmetric partner 
in connection to the Higgs physics given the large top Yukawa coupling, the bottom squark 
(sbottom b) sector is also of great interest. The left-handed sbottom mass is related to 
the left-handed stop mass since they are controlled by the same soft SUSY breaking mass 
parameter [4, 5]. In the region of large ratio of the Higgs vacuum expectation values: 
tan/3 = V 2 /V 1 , the bottom Yukawa coupling is large and there could be large corrections 
to the Higgs physics from the sbottom sector as well [6]. Although the LHC program has 
been carrying out a rather broad and impressive SUSY search plan, many searches are still 
under strong assumptions for the sake of simplicity. The current sbottom search mainly 
focuses on the direct decay to the lightest supersymmetric particle (LSP) b —)■ bxi, with 
bb -|- being the dominant search channel. With the data collected at the LHC 7 and 8 
TeV, a sbottom with mass up to 700 GeV has been excluded in this channel [7-9]. Even in 
the parameter space with highly degenerate sbottom and LSP masses [10-14], a sbottom is 
excluded with mass up to about 255 GeV [9]. Other decay channels including the cascade 
decay via the next-lightest supersymmetric particle (NLSP) b —?■ fey® ^hxi, bZxi [15, 16] 
and b —?■ tWxi [16-19] have also been considered with a 100% branching fraction each, 
with considerably weaker limits. 

In realistic MSSM scenarios, there are often more than one significant decay modes to 
be present. Two prominent examples stand out; A left-handed sbottom in the Wino-NLSP 
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scenario may have a decay b —)■ 6 x 2 with branching fraction as high as 30% — 40%, along 
with the leading decay b tyf; Similarly, a right-handed sbottom in the Higgsino-NLSP 
scenario may have the leading decay mode b bxi with branching fraction only 40%—60%, 
along with a sub-leading decay b —)■ tx^ of 20% — 30%. Those additional channels dilute 
the leading signals currently being searched for at the LHC, and signihcantly weaken the 
sbottom search limits when assuming 100 % branching fraction for a given search channel. 
On the other hand, the combination of the multiple decay modes offers alternative discovery 
channels for sbottom searches, that must be properly taken into account. 

In this paper, we present the sbottom decays in a few representative SUSY mass sce¬ 
narios. We then analyze the sbottom pair production signal with one sbottom decaying 
via b —)■ bxi, &X 2 ) the other one decaying via b —)■ txf- With the subsequent decay of 
X 2 ^ Zxi, hxi and Xi ^ ^^Xi, we study the reach of those signals at the 14 TeV LHC 
with 300 fb”^ integrated luminosity. Because of the similarity of the final state signatures 
and potential correlation of the left-handed soft mass, the sbottom and stop signals are 
combined to obtain the optimal sensitivity for the same SUSY parameter region. We find 
for a left-handed bottom squark, a mass up to 920 GeV can be discovered at 5a significance 
for 250 GeV < m^o < 350 GeV, or excluded up to 1050 GeV at the 95% confidence level 
for the h channel (p > 0); similarly, the bottom squark can be discovered up to 840 GeV, 
or excluded up to 900 GeV at the 95% conhdence level for the Z channel (/i < 0), the top 
squark reach is close to that of the bottom squark. The sbottom and stop signals in the 
same SUSY parameter scenario are combined to obtain the optimal sensitivity, which is 
about 150 GeV better than the individual reach of the sbottom or stop. For a right-handed 
bottom squark with the channel bb* —)■ bxi, tx^ —t tbW + we find that a mass up 
to 880 GeV can be discovered at 5a significance, or excluded up to 1060 GeV at the 95% 
conhdence level. 

The rest of the paper is organized as follows. In Sec. H, we briehy present the sbottom 
sector in the MSSM and introduce the mass and mixing parameters. We then calculate 
the sbottom decays for various neutralino/chargino mass spectra. Assuming one decay 
channel dominant at a time, we summarize the current LHG stop and sbottom search 
results from both ATLAS and GMS experiments. In Sec. HI, we investigate the reach 
of the sbottom signal with mixed decay channels at the 14 TeV LHG. We combine the 
left-handed sbottom and stop signals for the same SUSY parameter region. In Sec. IV, we 
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summarize our results. 


II. MSSM SBOTTOM SECTOR 

We work in the MSSM and focus primarily on the third generation squark sector. We de¬ 
couple other SUSY particles: the gluino, sleptons, and the hrst two generations of squarks. 
We also decouple the non-SM Higgs particles by setting Ma large. Besides the third gener¬ 
ation squarks, the other relevant SUSY states are a Bino (with a soft SUSY breaking mass 
Ml), Winos (with a soft SUSY breaking mass M 2 ), and Higgsinos (with bilinear Higgs 
mass parameter /i). Up on the mass diagonalization, they form neutralinos (X 1234 ) 
charginos {xt 2 )- 

A. The sbottom sector 

The gauge eigenstates for the the third generation squark sector are ii, bL, in, b^, where 
the left-handed states form a SU(2 )l doublet with the soft SUSY breaking mass M^sq, and 
the right-handed states are SU(2)i singlets with soft SUSY breaking masses M^su, M^sd- 
For the sbottom sector, the mass matrix in the basis of {b^, bn) is [4, 5] 

MisQ + ml + nibAb 

rribAb Misn + ml + 

where 

A 4 = (“^ + ^ sin^ 9w) cos 2/3M|, Aj^ = ^ sin^ 9w cos 2/3M| ( 2 ) 

are the contributions from the SU(2 )l and U(l)y D-term quartic interactions. The trilinear 
soft SUSY breaking coupling Ab leads to the off-diagonal term Ab = Ab — /itan/3, that 
induces the mixing between left-handed and right-handed sbottom states. The lighter and 
heavier mass eigenvalues will be denoted as respectively. 

The left-handed mass parameter M^sq also controls the mass of the lighter stop. Since 
the stop sector provides the dominant contribution to the Higgs mass corrections, “natu¬ 
ralness” argument prefers a relatively lower value of the stop mass. It is thus reasonable 
to consider rubAb, s^^icl the lighter sbottom mass eigenstate is mostly left- 

handed bi ~ The left-handed sbottom couples to a bottom quark and a neutralino (or a 
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top quark and a chargino) mainly through SU(2) l gauge coupling and top Yukawa coupling, 
depending on the components of the neutralinos/charginos. Although the sbottom correc¬ 
tions to the Higgs mass are small compared to that from the stop, there can be signihcant 
modihcation to the Higgs couplings, especially the bottom Yukawa coupling [6, 20]. 

A right-handed sbottom fe/j, on the other hand, couples to the U(l)y gaugino and 
Higgsinos only via the hyper-charge and Yukawa coupling. Its mass is determined by 
M^sd- consider the situation when it is light. 

B. Sbottom decays 

The most commonly studied channel in experimental searches is the case bi —?■ bxi 
with a branching fraction of 100%. This is true for the case with the Bino-LSP and the 
sbottom-NLSP, or the case with the stop-NLSP but mi < m£ + Mw, or the case with the 
Wino-NLSP for a right-handed sbottom. In a more general ground, sbottom decays lead 
to a much richer pattern. 

1. The decay of bL 

For a more general electroweakino spectrum, other decay channels may appear or even 
dominate, as analyzed in detail in Ref. [21]. We hrst consider the case of the lighter 
sbottom being mainly left-handed bi ~ bi- The mass spectrum of sbottom and gaugino 
would influence severely the decay modes of sbottom. We discuss the (mainly left-handed) 
sbottom decay in details in the two general situations with a Bino-LSP: 


mi^ > Ms > Ml 

(Wino-NLSP), 

(3) 

m-bi > Ihl > Ml 

(Higgsino — NLSP). 

(4) 


The more involved cases when both Winos and Higgsinos are below the sbottom mass 
threshold 


mi^ > IpI > Ms > Ml (Wino - NLSP/Higgsino - NNLSP), (5) 

mi^ > Ms > IpI > Ml (Higgsino - NLSP/Wino - NNLSP), (6) 

are also included when distinct features are present (sometimes referred as mixed NLSP’s). 
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A'h = 300 GeV, Ml = 150 GeV /i = 300 GeV, Mi = 150 GeV 



fj, = 450 GeV, Ah = 300 GeV, Mi = 150 GeV Ah = 450 GeV, = 300 GeV, Ah = 150 GeV 



FIG. 1: Branch fractions of the left-handed sbottom decay versus its mass in four different 
cases: (a) > M 2 > Mi, Wino-NLSP; (b) > |/i| > Mi, Higgsino-NLSP; (c) > |//| > 

Ma > Ml, Wino-NLSP/Higgsino-NNLSP, and (d) > M 2 > |/i| > Mi, Higgsino-NLSP/Wino- 

NNLSP. Here we have adopted tan/3 = 10. 

We illustrate the sbottom decay in Fig. 1 for these four different situations. Each corre¬ 
sponds to a different mass spectrum of gaugino and sbottom for a Bino-LSP. The usually 
considered channel bxi is suppressed, if other channels are open, since the bino U(l)y 
coupling is smaller than the wino SU(2)i coupling or top Yukawa coupling. In Fig. 1(a), 
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A'h = 300 GeV, Ml = 150 GeV /i = 300 GeV, Mi = 150 GeV 



FIG. 2: Branch fractions of the right-handed sbottom decay versus its mass for (a) > M 2 > 

Ml, Wino-NLSP and (b) > \^\ > Mi, Higgsino-NLSP. Here we have adopted tan/? = 10. 

where hi txf and hi —)■ 6 x 2 open while the Higgsinos-like neutralinos/charginos are 

decoupling (|yu| > Mi^ > M 2 > Mi), sbottom decays dominantly into txf and hy^. Con- 
trarily, in Fig. 1(b), we decouple Wino-like gaugino while leaving the channel containing 
Higgsino-like gaugino opening (M 2 > M^^ > |/i| > Mi), hi txf will soon dominant 
over other possible channels when the phase space is open due to the large top Yukawa 
coupling, hi —)■ 6 x 2,3 suppressed due to the relatively small bottom Yukawa coupling. 
Here we have adopted tan/? = 10. For a larger value of tan/?, 6 x 256 x 3 channels will be 
relatively more important. For more complicated situation, in the lower two panels, we 
consider the cases of M^^ > |/i| > M 2 > Mi (Fig. 1(c)) and M^^ > M 2 > |/i| > Mi 
(Fig. 1(d)). In both cases, sbottom decays dominantly into Higgsino-like chargino, then 
Wino-like chargino and at last Wino-like neutralino. Other channels are highly suppressed 
since the U(l)y coupling and bottom Yukawa coupling are much smaller. 

A special remark is in order. Although 6 l and II share the same soft mass parameter 
M 3 SQ, the large mixing between ii — in due to the large trilinear soft SUSY breaking At 
often drags the mass of the (mixed) stop below that of the (mainly left-handed) sbottom. 
The decay hi ^ Wti usually dominates once it is kinematically open. However, the about 
decay patterns still hold as long as M^^ < M^^ -|- mw. 
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2. The decay of bji 


For the bn, the usually considered channel bxi is the dominant mode. We present the 
branching fractions of b^ in Fig. 2, for (a) the Wino-NLSP and (b) the Higgsino-NLSP. 
We see that the channel bi —?■ bxi in the Wino-NLSP scenario is almost 100%, since 
the right-handed squark has no SU(2)i coupling. However, this channel in the Higgsino- 
NLSP scenario presents a branching fraction about 40% — 60%, followed by the channel of 
bi —t txf about 20% — 30%, due to the coupling effects of the right-handed squark to Bino 
or Higgsino is U(l)y and bottom Yukawa, respectively. 

C. Current bounds from LHC 

Searches for direct stop and sbottom pair production have been performed at both 
ATLAS and CMS, with about 5 fb“^ data at s/s = 7 TeV and about 20 fb~^ data at 
^/s = 8 TeV [7-9, 15-19, 22-32]. The current reach for the stop is slightly worse than 
that of the sbottom, which has been summarized in Ref. [21]. The current searches for 
the sbottom mainly focus on the decay channel bi —?■ bxi assuming 100% decay, and the 
sbottom mass up to 620 (700) GeV is excluded at 95% C.L. for a massless LSP with two 
b plus hnal states based on ATLAS (CMS) analyses [7, 8]. For small mass splitting 
between sbottom and the LSP: ~ nib, monojet plus search excludes a sbottom 

mass up to about 255 GeV [9]. 

Sbottom searches for b —)■ bx 2 , X 2 Xi^ with 100% decay branching fraction have been 
performed at ATLAS [15] and the null search results exclude the sbottom masses between 
340 and 600 GeV for m^o = 300 GeV and m^o = 60 GeV. For b -x- bx 2 , X 2 ^ Xi^ 

100% decay branching fraction, GMS searches exclude sbottom masses up to 450 GeV for 
LSP masses between 100 to 125 GeV and m^o — m^o = 110 GeV [16]. Sbottom searches 
for b —)■ txf, xf Wxi with 100% decay branching fraction have been performed at 
both ATLAS and GMS [17, 18]. The sbottom mass limit by ATLAS is about 440 GeV 
for m^± < — rrit [17], and the GMS limits for those channels are about 50 to 100 GeV 
stronger [16, 18, 19]. We summarize the current search bounds in Table 1. 
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Decay channels 

Mass bounds mj^ 

BR 

Assumptions 

bi (ATLAS [8]) 

620 GeV 

520 GeV 

100% 

60% 

m^o <120 GeV 

m^o < 150 GeV 

bi (ATLAS [9]) 

255 GeV 

100% 

mi — m^o ~ rufe 

h ^ (CMS [7]) 

700 GeV 

100% 

Small m^o 

bi bxl hhx\ (ATLAS [15]) 

340 - 600 GeV 

100% 

m^o = 300 GeV 

m^o = 60 GeV 

Xi 

6i ^ bxl ^ bZxl (CMS [16]) 

450 GeV 

100% 

100 GeV < m^o <125 GeV 

A.1 

m^o — m^o = 110 GeV 

2 C 2 

h txi (ATLAS [17]) 

440 GeV 

100% 

m^± <mi — mt 


575 GeV 

100% 

150 GeV < m + < 375 GeV 

Xi 

m o = 50 GeV 

Al 

bi ^ txT (CMS [16]) 

575 GeV 

100% 

25 GeV < m^o < 150 GeV 

Al 

m 0 

^1 —05 

m ± 


525 GeV 

100% 

25 GeV < m^o < 200 GeV 

Al 

""I = 0.8 

m + 

Al 

6i ^ txT (CMS [18]) 

500 GeV 

100% 

m 0 

= 0-5 (0.8) 

Al 

6i ^ txi (CMS [19]) 

550 GeV 

100% 

m^o = 50 GeV 

Al 


TABLE I: Current mass bounds on the sbottom from the direct searches at the LHC. 


III. LHC ANALYSIS 

In this section, we stndy the collider phenomenology of the light sbottom at the 14 
TeV LHC. The key point in this paper is to explore the mixed decay channels according 
to the mass hierarchies beyond the common assnmption of 100% branching fraction of a 
given channel. Including those channels listed in Table I with realistic branching fractions 
would help increase the overall sensitivity, but we did not repeat the analyses. We note 
that Ref. [33] also exploited the mixed decays to search for stop. They introduced a new 
variable “topness” for the top-rich signal events to help efficiently reduce the top pair 
backgrounds. 
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Ml 

M 2 

M 3 SQ 

At 


tan /3 

m^o 

Xi 

m^o 

A2 

m ± 
Xi 



rrih 

BPl 

150 

300 

650 

2950 

-h2000 

10 

152 

320 

320 

640 

650 

125 

BP2 

150 

300 

650 

2950 

-1300 

10 

150 

320 

320 

640 

630 

125 


TABLE II: MSSM parameters and mass spectrum of SUSY particles for the two benchmark 
points. All masses are in units of GeV. 

A. Signature of bi hi 

We consider the scenario with the low energy mass spectrum containing a light sbottom 
(mostly left-handed), a Bino-like LSP and Wino-like NLSPs, as shown in Fig. 1(a). Two 
typical benchmark points for both signs of p are listed in Table. II. Other soft SUSY 
breaking parameters are decoupled to be 2 TeV, and At is set to be large such that the 
SM-like Higgs is around 125 GeV. The value of /i is chosen such that X 2 dominantly decays 
to hxi for p > 0 and to Zxi for /i < 0^. The decay channels and the corresponding 
decay branching fractions for 61 , U, as well as X 2 xt listed in Table. III. The 
conventional channel bi —)■ bxi is highly suppressed, with only about 2 % branching fraction, 
which dramatically weakens the current experimental search limit. The decay channels of 
bi —^ bx 2 and bi —)■ txi are comparable and dominant instead. In particular, with one 
sbottom decaying to X 2 and one sbottom decaying to yf, bibl pair production leads to 
interesting hnal states of bbWW+h/Z+^rp. Note that unmixed decays of bib^ bbhh+^j,^ 
bbZZ + ttWW + $rp have been studied at the LHC [15-19], assuming 100% decay 
branching fractions. Given the more realistic branching fractions of about 40% for bi —?■ bx% 
and about 60% for bi —)■ txi ^ the collider limits for those channels will be relaxed. Including 
all the mixed and unmixed channels can further increase the collider reach for the sbottom. 

The stop decay has been studied in detail in Ref. [21]. For the two benchmark points 
listed in Table. H, the conventional decay channel H txi is highly suppressed. U —?■ bxi 
is dominant with branching fraction of about 70%. ti tx2 is subdominant with a 
branching fraction of about 27%. With one stop decaying to X2 one stop decaying to 
Xi, iiil pair production provides the same hnal states as the sbottom case. 

^ Note that X 2 ^Xi is not always dominated for ^ < 0, as pointed out in Refs. [34, 35]. We have chosen 
the value of ^ in the ^ < 0 case to guarantee the Z channel dominance. 
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Decay Channel 

BR 

Decay Channel 

BR 

Decay Channel 

BR 

BPl (p > 0) 

bi bxi 

2% 

ii txi 

2% 

X 2 hxl 

97% 

bi bx2 

39% 

ii tx2 

27% 

^ Zxl 

3% 

h txi 

59% 

ii bxt 

71% 

xt W^Xi 

100% 

BP2 (p < 0) 

bi bxi 

2% 

ii txi 

2% 

X 2 hxl 

6% 

bi bx2 

39% 

ii tx2 

27 % 

X^ ^ ^X? 

94% 

hi txi 

59% 

ii bxt 

71% 

xt ^ w^x^i 

100% 


TABLE III: Decay channels and the corresponding branching fractions of bi, ii, X 2 xt 
the two benchmark points, which correspond to the cases of /r > 0 and /r < 0. 

The two benchmark points listed in Table. II are only for illustration whenever instruc¬ 
tive. In our following analyses, we perform a broad scan over the mass parameter space. 

• from 400 to 1075 GeV with a step size of 25 GeV, corresponding to from 

about 350 GeV to about 1085 GeV and from about 367 GeV to about 1090 GeV. 

• Ml is scanned from 3 GeV to 700 GeV, in the step of 25 GeV. 

• M 2 is hxed to be M 2 = Mi -|- 150 GeV. 

• We further require > m^± -|- rrit such that hi —)■ txi can be open. 

In our phenomenological studies, we dehne the basic observable objects as 

• Jet: 

\7],\ < 2.5, 4 > 25 GeV, > 0.8. (7) 

where is azimuthal angle between the jet and missing transverse energy. 

• Lepton: 

\t]i\ < 2.5, > 20 GeV, Ai?^^ > 0.4. (8) 

Where the Ai?^^ is the distance in the 0-?7 plane: Ai? = A0^ -|- Arf, between the 

lepton and the jet satisfying Eq. (7). 

To be as realistic as possible, both the signal and the background samples are generated by 
MadGraph 5 [36], passed through Pythia 6 [37] for the fragmentation and hadronization. 
We further perform the detector simulation through Delphes 3 [38] with Snowmass Delphes 
No-Pile-up detector cards [39]. 
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1. The Case of > 0; final states with a Higgs 


In the case of /i > 0, the leading signal nnder consideration for the pair prodnction of 
sbottom, with hi 6 x 2 bhxi and 6 * —)■ txf —t bW^ W~Xi, is 

bib*i bb WW bbbb jj$T- 

The signal contains fonr 6 -jets, two light flavor jets, one isolated lepton (e or p), and 
large missing energy. The stndy of the same final state from stop decay can be fonnd 
in Ref. [21]. The dominant backgrounds will be from ff-Tjets and tibh with large cross 
sections and similar final states. While tih is an irreducible background, the production 
cross section is relatively small. Other SM backgrounds include tfW, tiZ and hhWW, with 
typically smaller cross sections. 

To select the signal of iii*i —)■ bb WW h $rp —)■ f bbbb jj we adopt the basic 
event selection 

• Nj > 4, > 40 GeV, W = 1- 

Beside these basic cuts, we further optimize the cuts and divide the events into signal 
regions on the following variables: 

• Missing energy which is the magnitude of the the missing transverse momentum, 
to be above 100, 120, 140, 160 180, and 200 GeV. 

• Hti the scalar sum of the jet transverse momentum of all surviving isolated jets: 
Ht = Xljets IPtI) to be above 400, 450, 500, 550, 600 GeV. 

• Mjr, the transverse mass, defined as the invariant mass of the lepton and missing 
energy: 

Pt(9) 
to be above 100, 120, 140, 160, 180, 200 GeV. 

• Nj, the multiplicity of all surviving isolated jets, being at least 4, 5 and 6 . 

• Nf,, the multiplicity of tagged 6 -jets, being at least 2, 3 and 4. 
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FIG. 3: Normalized distributions of (left panel) and Ht (right panel) for the signal bib'l (red 
curves), (blue curves) —^ bbWWh$j’ —)• i bbbb jj $rp after basic cuts with = 637 GeV, 
= 646 GeV, as well as SM backgrounds at the 14 TeV LHG. 

The normalized distributions of ^ind Ht are shown in Fig. 3. As expected, the signal 
process has larger $t from the missing neutralino-LSP than the background processes, 
which is typically bounded by mw/‘^ due to the primary contribution W Given 

the relatively large sbottom mass, the signal process typically has larger Ht than the SM 
backgrounds as well. 

In Table. IV, we list the cumulative cut efficiencies after different levels of cuts, as well 
as cross sections before and after cuts for both the sbottom and stop signals as well as the 
SM backgrounds for the benchmark point listed in Table II for p > 0. The cross section 
for each process is normalized to their theoretical values including NLO QCD corrections 
[40-46]. The background processes are significantly suppressed after strong $t-i Ht, Mt 
cuts. The leading background left is ti, followed by tibb. We scan over the combinations 
of the signal regions, to select the optimal combination which gives the best significance 
for each mass grid point, including 10% systematic uncertainty. At -y/s = 14 TeV with 
300 fb“^ integrated luminosity, the signihcance could reach about 17a (14cr) for bi (G) of 
about 640 GeV. 

Signal signihcance contours are shown in Fig. 4 with the 5a discovery reach (black 
curve) and 95% G.L. exclusion limit (red curve) for 14 TeV LHG with 300 fb~^ integrated 
luminosity. Fig. 4 (a) shows the — m^o plane. We hnd that 5cr discovery can reach 
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Process 

a (fb) 

Basic 

$T > 

Ht > 

Mx > 

Nj> 

Nb> 

a (fb) 



cuts 

200 GeV 

500 GeV 

160 GeV 

5 

2 

after cuts 

hbi 

13 

39% 

17% 

14% 

5.8% 

4.3% 

2.7% 

3.4 X 10-1 

iih 

10 

39% 

18% 

16% 

5.9% 

4.4% 

2.9% 

2.9 X 10-1 

tt 

260,000 

14% 

0.24% 

7.4 X lO”"^ 

1.7 X 10“® 

9.3 X lO-’^ 

2.4 X 10-"^ 

6.3 X 10-2 

tibb 

2,300 

24% 

0.6% 

0.3% 

3.5 X 10-5 

2.3 X 10-5 

1.2 X 10-5 

2.8 X 10-2 

tih 

100 

31% 

1.2% 

0.8% 

5.8 X 10-5 

3.4 X 10-5 

1.9 X 10-5 

2.0 X 10-3 

tiZ 

230 

30% 

1.2% 

0.8% 

6.6 X 10-5 

3.9 X 10-5 

9.8 X 10-5 

2.2 X 10-3 

tiW^ 

224 

25% 

1.2% 

0.7% 

4.8 X 10-5 

2.3 X 10-5 

6.3 X 10-5 

1.4 X 10-3 




vG = 14 TeV j Ldt = 300 fb”^ 





s 

x/n+(io%s 

- - 17 (14) for h {ii) 




TABLE IV: Cut efficiencies and cross sections before and after cuts for the signal —)■ 

bbWWh$T ^ bbbb jj for BPl listed in Table II for fi > 0, as well as SM backgrounds at 
the 14 TeV LHC. The significance is obtained for f Ldt = 300 fb~^ with 10% systematic error 
combining both sbottom and stop signals. 

about 750 GeV for hi when Xi is almost massless and reach about 920 GeV when Xi is 
about 200 GeV to 300 GeV. The 95% G.L. exclusion reach is about 100 GeV better. The 
reach for the stop with the same hnal states can be found in Ref. [21], with results being 
very similar. 

Since the (mostly left-handed) sbottom and stop have the same undistinguishable final 
states with their masses controlled by the same parameter M^sq, we present the combined 
reach of stop and sbottom in Fig. 4 (b) in M^sq ~ The 5a discovery reach 

in Mssq increases to be 820 GeV for a massless LSP, and 1080 GeV for m^o ~ 300 GeV. 
The masses up to 980 GeV can be excluded for a massless LSP, and the masses up to 1180 
GeV can be excluded for m^o ~ 300 GeV at 95% G.L. 

We would like to reiterate that the mixing in sbottom and stop sectors governs the 
mass spectrum of the sbottom and stop. Small mixing in the sbottom sector is always a 

^ The mass difference between the stop and sbottom does not affect the combination of the stop and 
sbottom signals, since the same cuts are used for both the stop and sbottom events. 
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FIG. 4: Signal significance contours for bib*,iii\ —)• bbWWh^x —I bbbb jj $j' final states for 
14 TeV LHC with J L dt = 300 fb“^ luminosity. The 5a discovery reach (black curves) and 
95% C.L exclusion limit (red curves) for the sbottom only are shown in the (a) — m^o plane, 

and in the (b) M^sq — plane for the combined reach for sbottom and stop. 


good approximation given the small bottom Yukawa coupling, while the mixing in the stop 
sector may be large enough to suppress the mass of the lighter stop further. In our cases 
(including the /i < 0 case discussed below), the right-handed stop is assumed to be very 
heavy (decoupled to be 2 TeV), which will result in a smaller mixing for a large range of 
At G [—4000,4000] GeV. Furthermore, even if a large mixing in stop sector gives a much 
lighter stop compared with the sbottom, this would potentially lead to a better signal 
in the stop sector. The combination of the stop and sbottom signals, however, does not 
depend on the mass difference between the stop and sbottom. In the parameter space that 
we are considering with relatively small stop and sbottom mass difference, both channels 
contribute signihcantly to the combined reach. In cases when the mass difference between 
the stop and sbottom is large, only one channel will contribute dominantly to the combined 
signihcance. 


2. The Case of fj, < 0." final states with a Z-boson 

For the case of /i < 0, the dominant decay channel of X 2 X 2 ^Xi instead [47]. The 
leading signal under consideration for the pair production of sbottom and stop with bi 
bx^ bZxi, K txf hW+ W-x\ and h txl ^ hW+Zx\, il bxf bW-Xi, is 
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then 


bibl,iiil bb WW Z ^ bb jjjj $t- 

The signal contains two 6-jets, four light flavor jets, two same flavor, opposite sign leptons, 
and large missing energy. The two leptons are used to reconstruct the Z boson, which will 
significant reduce the SM backgrounds. The dominant background is ti plus one or two 
additional QCD jets. 

We impose the basic event selection cuts as the previous case. We again optimize the 
cuts and divide the events into signal regions: 

• $j’ to be above 100, 120, 140, 160 180, and 200 GeV. 

• Ht to be above 400, 450, 500, 550, 600 GeV. 

• Mt 2 , the lepton-bashed transverse mass [48-50]: 

M^2(PT^PT^PT"^)= . min . {max{M^(p^^p”f),MT(p^^p?^r)}} (10) 

--.miss _i__miss _.miss ’ ’ 

1 “rPj- 2 —Pj- 

to be above 75, 80, 85, 90 GeV. 

• AMi^ = \Mu — Taz\., being less than 10 GeV. 

• Nj being at least 4, 5 and 6. 

• Nb being at least 1 to suppress the enormous QGD backgrounds with light falvor 
jets. 

The normalized distributions of Mt 2 for both the sbottom and stop signal, as 

well as the SM backgrounds are presented in Fig. 5. The distributions for the signal 
typically extend to larger values. The Mt 2 {Pt ^ Pt > Pt'^^) distributions for SM backgrounds 
with the lepton pair coming from leptonic W decay are cut off at mw, while the signal 
as well as bbZZ background have much flatter Mt 2 distributions. Note that while the 
distribution of bbZZ background is similar to that of the signal, the overall cross section 
for bbZZ is negligibly small. 

Another interesting variable for the sbottom case is which is related to if the b 
jet and the lepton pair from the same sbottom cascade decay chain bi bx 2 bZxi can 
be identihed. While we will not use it for event selection in our analyses. Mm distribution 
could provide information on as well as m^o if a sbottom signal is discovered. 
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FIG. 5: Normalized distributions (left panel) and Mt 2 (right panel) for the signal bib* (red 
curve), iiil (blue curves) — bbWWZ$rp l'^i~ bb jjjj $rp after basic cuts with = 637 
GeV, = 634 GeV, as well as SM backgrounds at the 14 TeV LHC. 

The advanced cuts and the corresponding cumulative cut efficiencies as well as the cross 
sections for sbottom and stop signal for BP2 with p < 0 and SM backgrounds before and 
after cuts are given in Table V. The dominant SM background is ti plus jets. A significance 
of about 12a {8.7a) can be reached for bi {ti) for the benchmark point at the 14 TeV LHC 
with 300 fb~^ luminosity, including 10% systematic error. 

Signal significance contours are shown in Fig. 6 with the 5cr discovery reach (black 
curve) and 95% C.L. exclusion limit (red curve) for 14 TeV LHC with 300 fb”^ integrated 
luminosity, in the (a) (b) mi^ ~'^xi and (c) M^sg — rn^o plane. For 

massless y?, sbottom (stop) masses up to 650 (680) GeV can be discovered and 720 (760) 
will be excluded at 95% C.L. if there is no signal over SM backgrounds being found. For 
moderate mass of Xi around 200 ~ 300 GeV, the 5a dicovery can reach up to 820 (840) 
GeV, and the 95% exclusion limit can go up to 890 (900) GeV for sbottom (stop). The 
combined reach of the stop and sbottom is shown in Fig. 6 (c) in M^sg versus m^o plane. 
About 980 GeV can be achieved in M^sg for the 5cr discovery reach and about 1025 GeV 
for the 95% C.L. exclusion. The experimental reach for the case of p < 0 is lower than 
that for the case of /i > 0. 
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Process 

a (fb) 

Basic 

$T > 

Ht > 

Mt2 ^ 

AMu < 

N,> 

Nb> 

a (fb) 



cuts 

175 GeV 

400 GeV 

90 GeV 

10 GeV 

4 

1 

after cuts 

hbl 

2.1 

32% 

17% 

16% 

5.8% 

5.3% 

5.3% 

4.2% 

8.8 X 10-2 

iJi 

1.8 

27% 

16% 

11.2% 

4.9% 

4.4% 

4.4% 

3.6% 

6.5 X 10-2 

tt 

33,000 

1.3% 

0.09% 

0.06% 

5.0 X 10“® 

4.9 X 10-^ 

4.9 X 10-^ 

3.7 X 10-^ 

1.2 X 10-2 

tiZ 

71 

11% 

0.25% 

0.16% 

5.8 X 10-5 

4.2 X 10-5 

4.2 X 10-5 

2.7 X 10-5 

1.9 X 10-5 

ttbb 

400 

3.2% 

0.20% 

0.12% 

1.4 X 10-5 

2.0 X 10-5 

2.0 X 10-5 

1.8 X 10-5 

6.9 X 10-4 

ttzz 

0.16 

16% 

0.86% 

0.64% 

0.31% 

0.27% 

0.27% 

0.18% 

3.0 X 10-4 

bbZZ 

2.3 

0.39% 

0.11% 

0.06% 

2.9 X 10-^ 

2.6 X 10-^ 

2.6 X 10-4 

2.1 X 10-4 

4.8 X 10-4 





= 14 TeV JL dt = 300 fb"i 







S — 1o 

^yB+{l0%B)^^ 

(8.7) for h (ii) 




TABLE V: Cut efficiencies and cross sections before and after cuts for the signal —>• 

bhWWZ^rp —>■ bb jjjj for BP2 in Table II for /i < 0, as well as dominant SM back¬ 
grounds at the 14 TeV LHC. The significance is obtained for J L dt = 300 fb“^ with 10% 
systematic error combining both sbottom and stop signals. 

B. Signature of bi r\j bR 

To complete our exploration for the sbottom signal, we consider another scenario with 
the low energy mass spectrum containing a light mostly right-handed sbottom, a Bino-like 
LSP and Higgsino-like NLSPs. Here, the sign of /i does not affect the decay modes of sbot¬ 
tom and neutralinos. The typical benchmark point is listed in Table. VI, the corresponding 
branching fractions are listed in Table. VII. Other soft SUSY breaking parameters are de¬ 
coupled by setting them to be at 2 TeV. In this scenario, the right-handed sbottom couples 
to the Bino and Higgsino through the U(l)y or the bottom Yukawa couplings, which re¬ 
sults in the sbottom dominantly decaying to bxi due to the large phase space, followed by 
the channel tx^ when it is kinematically open. We will focus on the signal reach of the 
sbottom pair production 

hbl bxi txf bxi tW^Xi ^ ^bb jj 

The SM backgrounds are somewhat similar to that of the fi > 0 case of left-handed sbottom 
with fewer jets. We also include vector bosons plus additional jets as another background 
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FIG. 6: Signal significance contours for hih*^iii\ —)• bbWWZ^j- — t^i~bbjjjj$j’ final states 
for 14 TeV LHC with f L dt = 300 fb~^ luminosity. The 5a discovery reach (black curves) and 
95% C.L. exclusion limit (red curves) are show in the (a) — m^o plane, in the (b) — m^o 

plane, and the combined reach in the (c) M^sq — 'm^o plane. 



Ml 

Ma 

MzSD 

At 

h 

tan (3 

m^o 

A.1 

^xt 


mh 

BPS 

150 

2000 

650 

2895 

300 

10 

145 

307 

635 

125 


TABLE VI: MSSM parameters and mass spectrum of SUSY particles for the the benchmark point 
in the case of Right-handed sbottom. All masses are in units of GeV. 

|51]. 

We scan over a broad mass parameter space: Mi from 3 GeV to 800 GeV in step 
of 30 GeV, M^sd from 400 GeV to 1180 GeV in step of 30 GeV, fi is fixed to be /r = 
Ml + 150 GeV. We further require that mi^ > m^± + nit so that the decay channel 
^xt is kinematically accessible. Since the hnal state particles are more stiff than the 
previous cases with cascade decays, we apply stronger basic cuts than before. The basic 
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Decay Channel 

BR 

Decay Channel 

BR 

BP3 

61 ^ 6 x 1 

58% 

xf ^ X?1V± 

100 % 

61 ^ txi 

18% 




TABLE VII: Decay channels and the corresponding branching fractions of bi and for the 
benchmark point, which corresponds to the case of Right-handed sbottom. 

event selection cuts are 

• Jet: 

\r]j\ < 2.5, pip > 40 GeV, > 0.8. (11) 

• Lepton: 

\pi\ < 2.5, P^rp > 30 GeV, > 0.4. (12) 

• at least three jets satisfying requirement Eq. (11), within which at least one fe-tagged, 
and exactly one lepton satisfying requirement Eq. (12). 

• the leading 6-jet px is required to be larger than 100 GeV since one 6-jet originates 
directly from a heavy sbottom decay. 

Besides the basic event selection cuts, we apply the same advanced event selection cuts 
in the signal regions (Lfr, $t-i Afr, Nj and Ni,) as in Sec. Ill A, and optimize them for 
different mass parameters. In Table VIII, we list the cross section before and after above 
cuts and also the efficiency after every cut for the benchmark point listed in Table VI. 

Signal signihcance contours are shown in Fig. 7 with the 5cr discovery reach (black 
curve) and 95% G.L. exclusion limit (red curve) for 14 TeV LHG with 300 fb~^ integrated 
luminosity, in the — m^o plane. For a large range of mass of Xi (from massless to about 
300 GeV), sbottom masses up to about 880 GeV can be discovered and 1050 GeV will be 
excluded at 95% G.L. if there is no further signal over SM backgrounds being found. The 
reach at lower mass of Xi is better than that of left-handed case, since lowering the mass of 
Xi will increase the pr of the 6-jet produced together with Xi, and this effect is suppressed 
in the left-handed case where the leading 6-jets is produced together with X 2 Xi which 
are always heavier than Xi- 
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Process 

a (fb) 

Basic 

Cuts 

t'T > 

200 GeV 

Ht > 

500 GeV 

Mj' > 

160 GeV 

N,> 

4 

Nb> 

1 

a (fb) 

after Cuts 

bibi 

9.7 

30% 

20% 

14% 

8.1% 

5.6% 

5.6% 

5.4 X 10-1 

ti 

260,000 

5.3% 

0.14% 

4.7 X 10"^ 

1.6 X 10"® 

8.1 X 10-^ 

8.1 X 10-^ 

2.1 X 10-1 

ttbb 

2,300 

13% 

0.4% 

0.2% 

3.7 X 10-5 

2.6 X 10-5 

2.6 X 10-5 

6.1 X 10-2 

tth 

100 

20% 

1% 

0.7% 

7.8 X 10-5 

5.2 X 10-5 

5.2 X 10-5 

5.3 X 10-5 

ttz 

230 

14% 

0.7% 

0.5% 

8.1 X 10-5 

4.5 X 10-5 

4.5 X 10-5 

1 X 10-2 

ttw^ 

224 

11% 

0.7% 

0.5% 

6.6 X 10-5 

3.4 X 10-5 

3.4 X 10-5 

7.6 X 10-5 

Vjj 

3.7 X 10^ 

4.8 X 10“^ 

2.9 X 10“® 

1.8 X 10"® 

2.9 X 10-^ 

1 X 10-^ 

1 X 10-^ 

3.8 X 10-2 



^/s = U TeV L 

= 300fb“^ 

S/^B + (10%H)2 


11.4 


TABLE VIII: Cut efficiencies and cross sections before and after cuts for the signal hih\ 
bhWW$rp —). £ 55 jj for BP3 in Table VI, as well as SM backgrounds at the 14 TeV LHC. 
The significance is obtained for J Ldt = 300 fb~^ with 10% systematic error. 



Mg [GeV] 


FIG. 7: Signal significance contours for 6i6| —)■ —>• bb jj final state for the 

right-handed sbottom in the TeV LHC with J L dt = 300 fb“^ luminosity. 

The 5(7 discovery reach (black curves) and 95% C.L. exclusion limit (red curves) are shown. 

IV. SUMMARY AND CONCLUSION 

In this paper, we stress the point that in a realistic situation in a generic MSSM, the 
sbottom decay can be far from 100% to a specific channel, as assumed in most of the 
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current studies and all the LHC sbottom searches, which is only true for the Bino-LSP and 
either the left-handed sbottom (or a nearly degenerate stop) being the NLSP, or the right- 
handed sbottom (or Wino) being the NLSP. On a more general ground, sbottom decays 
lead to a much richer pattern. The inclusion of the other decay channels will significantly 
weaken the current sbottom search limits and in the mean time open new decay modes for 
alternative discovery channels for sbottom searches. 

We studied in detail the sbottom decay patterns in a few representative SUSY mass 
scenarios. For the left-handed sbottom, we found that 

(1) in the Wino-NLSP case, see Fig. 1(a), BR(fei —> bxl) ~ BR(fei —)■ txf ) ~ 50% while 
BR(6i ^ bx°i) ~ 2%. 

(2) in the Higgsino-NLSP case, see Fig. 1(b), bi —)■ dominates while fei —)■ bxi 2 3 

all suppressed. 

(3) in mixed NLSP cases, see Figs. 1(c) and 1(d), BR(6i ^ bx^) ~ BR(6i ^ txf) ~ 
BR(fei txf) ~ 30% and BR(fei —)■ bxf) ~ 3% when |/i| > Mx, while BR(fei 
txf) ~ BR(6i txf) ~ 30% and BR(fei —)■ bxf) < 10% when M 2 > \fi\. 

For the right-handed sbottom, see Fig. 2, decays of bi —)■ bxi dominates for the case of 
Bino-LSP with Wino-NLSP. In the case of Bino-LSP with Higgsino-NLSP, however, the 
branching fraction of bi bxi is reduced to about 40%—60%, while bi txf is about 
20—30%, followed by bi —)■ bxf^ of about 10% each. 

We analyzed in detail the sbottom pair production signals with the mixed decay chan¬ 
nels. We focus on the search sensitivity at the 14 TeV LHC with a 300 fb“^ integrated 
luminosity. We scanned over a large SUSY mass parameter region and performed semi- 
realisc detector simulations. For the left-handed sbottom bi pair production, we focused 
on the scenario of Bino-LSP with Wino-NLSP. With one sbottom decaying via b —)■ bx 2 
and the other sbottom decaying via b —)■ txf, we found that 

• With X 2 ^X 2 if' > 0) xf W^Xi, fhe leading signal is the bbbb jj I -f- 
hnal state. From Fig. 4(a), we see that a 5cr discovery can be made up to 920 GeV, 
and the 95% C.L exclusion limit can reach up to 1050 GeV for this Higgs channel. 
The reach of the combined sbottom and stop signals of the same final states is about 
120 GeV higher, as shown in Fig. 4(b). 
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• With X2 ^ Zx 2 (h < 0) and xt ^ W^Xi, we studied the reach of bb jjjj ££ + $t 
hnal state. As seen from Fig. 6, a 5cr discovery can be made up to 840 GeV, and 
the 95% C.L exclusion limit can reach up to 900 GeV for the Z channel. The 5a 
discovery potential of the combined sbottom and stop signals can reach up to 980 
GeV, and the 95% exclusion limit is about 1025 GeV. 

We also studied the signal for the right-handed sbottom bji in the scenario of Bino-LSP with 
Higgsino-LSP. With one sbottom decaying via b —)■ bxi, and the other sbottom decaying 
via b —)■ txf, we found that the reach of bbjj£ + hnal states can lead to a 5a discovery 
up to 900 GeV, and the 95% C.L exclusion limit up to 1060 GeV, as shown in Fig. 7. 
Including the other commonly studied channels, bxibxi, ^^^1 txiixt would help 

increase the overall search sensitivity, but we did not repeat the analyses as listed in Table 
I. 
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